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We obtain hydrodynamic equations describing a fluid consisting of chiral molecules or a suspension 
of chiral particles in a Newtonian fluid. The stresses arising in a flowing chiral liquid have a 
component forbidden by symmetry in a Newtonian liquid. For example, a chiral liquid in a Poiseuille 
flow between parallel plates exerts forces on the plates, which are perpendicular to the flow. A generic 
flow results in spatial separation of particles of different chirality. Thus even a racemic suspension 
will exhibit chiral properties in a generic flow. A suspension of particles of random shape in a 
Newtonian liquid is described by equations which are similar to those describing a racemic mixture 
of chiral particles in a liquid. 



Equations of hydrodynamics express conservation of 
mass, momentum and energy, and can be written as 

dtp + diJi^O, (la) 
^tP^ + d,U,j ^ 0, (lb) 
dtE + d^JE^O. (Ic) 

Here dt and di denote time and spatial derivatives, p, P, 
and E are correspondingly the densities of mass, momen- 
tum and energy, and J, J_b and 11 are the flux densities of 
mass, energy and momentum (we indicate vector quan- 
tities by bold face symbols and second rank tensors by 
hats). The flux densities can be expressed in terms of 
the hydrodynamic variables: the pressure p(r,t), tem- 
perature T(r,t) and the hydrodynamic velocity v(r,i), 
which we define via the equation 

pv = J = P. (2) 

To lowest order in spatial derivatives we have fT] 

n.y = pv^Vj +pSij - r/V^j - C6ij divv, (3) 

where Vij = djVi + diVj — ^5ij divv is the rate of shear 
strain, and rj and ^ are the first and the second viscosities. 
This leads to the Navier-Stokes equations, which should 
be supplemented by the equation of state of the fluid 
and the expression for the energy current in terms of the 
hydrodynamic variables. 

For a dilute suspension of particles in a Newtonian liq- 
uid, the basic hydrodynamic equations need to be sup- 
plemented [I] by the conservation law for the current of 
suspended particles, 

dtn + v ■ Vn + divj = 0. (4) 

Here n{r,t) is the density of suspended particles, and 
j(r,i) their flux density (relative to the fluid). To linear 
order in the gradients of concentration, temperature and 
pressure the latter can be written as 

j = -DVn - uXtVT - nXpVp. (5) 



Equation (|2| remains unchanged and can be considered 
as a definition of the hydrodynamic velocity v, which 
is, generally speaking, different from the local velocity 
u(r, t) near an individual particle of the suspension. 

There are corrections to the flux densities of various 
quantities, which are higher orders in spatial derivatives 
of the hydrodynamic variables (for a review see, for ex- 
ample, Refs. Moreover, there are nonlocal correc- 
tions to the Navier-Stokes equations, which can not be 
expressed in terms of higher order spatial derivatives of 
hydrodynamic variables [H O El |7j. 

Several studies focused on the effects of chirality on the 
motion of suspended particles in hydrodynamic flows [5J 
H m HU [H III m HSl US] . it was shown that non-chiral 
magnetic colloidal particles can self-assemble into chiral 
colloidal clusters [T7] . 

In this article, we develop a hydrodynamic description 
for the case of a suspension containing both right-handed 
and left-handed chiral particles in a centrosymmetric liq- 
uid. We show that in this case the corrections to the 
Navier-Stokes equations contain new terms, which are 
associated with the chirality of the particles. The sig- 
niflcance of these corrections is that they describe new 
effects, which are absent in the case of centrosymmetric 
liquid. Since certain types of hydrodynamic flow lead 
to separation of particles with different chirality, these 
corrections are important even in initially racemic sus- 
pensions of chiral particles. For simplicity we consider 
the case of incompressible fluids. 

In a given flow an individual particle of the suspension 
undergoes a complicated motion which depends on the 
initial position and orientation of the particle. The hy- 
drodynamic equations can be written for quantities which 
are averaged over the characteristic spatial and temporal 
scales of such motion. 

In the presence of chirality the following contribution 
to momentum flux density is allowed by symmetry: 

H^^ = n''''ar][d,ujj+djuj,] + r]ai[uj,djn''''+ujjd,n''% (6) 

where (x'i(r) — ^ CijkdjVki'r) is the flow vorticity, and 
n'^^ = (n+ — n_) is the chiral density, with n+ and 



n_ being the volume densities of right- and left-handed 
particles respectively. Equations (l][6) should be supple- 
mented by the expression for the chiral current, defined 
as the difference between the currents of right- and left- 
handed particles. Separating it into the convective part, 
vn^'* , and the current relative to the fluid, j'^'^, we write 
the continuity equation as 



dtUch + div(vn'='') -I- divj ''^ ^ 



(7) 



Besides the conventional contribution given by Eq. ([5| 
with n replaced by n'^^, the chiral current contains a con- 
tribution, j which depends on the flow vorticity: 



(8) 

where n = + n_ . The contributions to jf'^ containing 
only nuji are not allowed as there should be no chiral 
current in rigidly rotating fluid. 

In Eqs. Q and (|8| we keep only the leading terms in 
the powers of diVj , or in the order of spatial derivatives of 
V. Although these terms are subleading in comparison to 
those in the conventional hydrodynamic approximation, 
they describe new effects which are absent in the latter. 

Note that according to the Navier-Stokes equations 

V^curl (v) = -{afcurl (v) + curl [(vV)v]}. (9) 
V 

Thus the first term in Eq. ^ arises either due to non- 
stationary or non-linear in v nature of the flow. In par- 
ticular, in stationary flows and to zeroth order in the 
Reynolds number V^curl(v) = and this term vanishes. 

In spatially inhomogeneous flows the suspended parti- 
cles rotate, generally speaking, relative to the surround- 
ing fluid. This gives rise to separation of particles of 
different chirality due to the propeller effect, and to the 
chiral contribution to the momentum flux, Eq. ([6]). 

The rotation of the particles relative to the fluid arises 
due to two effects: 

i) In the presence of the spatial dependence of vorticity, 
Wi(r), the angular velocity of a particle is different from 
Li)i{r). This results in Eq. ^ and the first term in Eq. ([8|. 

ii) A non-uniform hydrodynamic flow induces orien- 
tational order in suspended particles similar to nematic 
order in liquid crystals. In the presence of flow vorticity 
orientation of particles induces their rotation with re- 
spect to the surrounding fluid. This contributes both to 
the chiral stress and the chiral flux. The latter contri- 
bution is described by the second term in Eq. (|8|. The 
contribution to the chiral part of the stress tensor asso- 
ciated with orientational order was discussed in Ref. |H1 

In most cases of practical importance the Reynolds 
number corresponding to the particle size R is small. In 
this regime the coefficients a, ai, (3, and /3i in Eqs. ([6| 
and ([s]) can be obtained by studying the particle motion 
in the surrounding fluid in the creeping flow approxima- 
tion |18l 119] . In this approximation the motion of parti- 
cle immersed in the liquid is of purely geometrical nature 



(see for example Ref. '50') . Dimensional analysis gives an 
estimate 



(10) 



where R is the characteristic size of the particles, and 
the dimensionless parameter x characterizes the degree 
of chirality in shape of the particles. 

The relative magnitudes of the different terms in 
Eqs. ^ and (|8| depend on the particle geometry. For 
example, particles with the symmetry of the isotropic he- 
licoid [IH] can not be oriented in a shear flow. Therefore 
the second term in Eq. ^ vanishes is this case. 

The degree of orientation of the particles can be ob- 
tained by balancing the characteristic directional relax- 
ation rate due to the Brownian rotary motion, ^ T/r/R^ 
with T being the temperature, with the rate of orienta- 
tion due to the shear flow, ^ Vij. Thus at small shear 
rates the degree of particle orientation is ^ Vij-qR^ /T . 
This leads to the estimate 



/3i ~ xvRVT. 



(11) 



The second term in Eq. ([8| is the leading term in the ex- 
pansion in the rotational Peclet number Pe ~ Vij-qR^ /T . 
At larger Peclet numbers terms of higher power in Vij 
should be taken into account. For Pc 1 the particle 
orientation becomes strong, and the corresponding con- 
tribution to the chiral current can be estimated as 



~'-ch 



J 



xRuj. 



(12) 



Equations ^ and ([s]) are written for the case when 
there is no external force acting on the particles, e.g. for 
a suspension of uncharged neutrally buoyant particles. 
In the presence of an external force F, there will be ad- 
ditional contributions to the chiral flux. The linear in 
F contributions can be constructed by contracting the 
antisymmetric tensor Cijk with the velocity Vi, force Fi 
and two derivatives di. For example, the following terms 
exist when F is constant: (F • V)u;, F x V^v, V(F • u}). 
These terms arise when the orientation of the suspended 
particles can be characterized by a polar vector. In this 
case the degree of particle orientation can be estimated 
as ~ RF/T. Thus the coefficients with which these terms 
enter the chiral current j '^^ are of order as xR^ /T- In the 
case when particles do not have a polar axis the degree 
of their orientation, and the corresponding contribution 
to the chiral flux are quadratic in F for small force. 

The chiral contribution to the stress tensor Eq. ^ 
leads to several new effects. Consider a Poiseuille flow 
of a chiral liquid between parallel planes separated by a 



distance d: Vx — —dxp{<P — 4y^)/8?7, 



(see 



Fig. [T]), with dxP being the pressure gradient along the 
flow. If the chiral density is uniform the chiral part of the 
stress tensor has only two non- vanishing elements, 11^^ = 
lizy = It describes a pair of opposing forces 
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FIG. 1: A chiral liquid in a Poiseuille flow between parallel 
plates exerts a pair of opposite forces on the plates = — F|, 
which are directed into ((g)) and out of (0) the page. A flow of 
a chiral liquid in a converging or diverging channel develops a 
helicoidal component of velocity directed into and out of the 
page, as shown at right. For a Newtonian fluid the flow lines 
(dotted lines) lie in the plane of the flgure. 



per unit area exerted by the liquid on the top and bottom 
planes. These forces are perpendicular to the flow, as 
shown in Fig. [l] Assuming n'^^R^ ^ 1 and using Eq. ( 10 ) 
the magnitude of the chiral force per unit area of the 
plane can be estimated as ~ xR'^xP- 

If rf'*^ is constant in space, then the volume force 
density generated by the chiral part of the stress ten- 



sor is /• 



ch 



arfsl oji. Then it is clear 



from Eq. (|9| that ff^ is generated only in nonstation- 
ary or nonlinear flows. In the special case of stationary 
Poiseuille flow the chiral part of the stress tensor does 
not generate a force density inside the fluid even at large 
Reynolds numbers. Thus the flow pattern is not affected 
by the fluid chirality. However, in a generic flow with 
converging or diverging flow lines the fluid chirality does 
affect the flow pattern. This is especially evident in flows, 
which have a mirror symmetry in the absence of chiral 
corrections. In these cases the chiral contribution to the 
stress tensor results in mirror asymmetric corrections to 
the flow velocity. For example a chiral liquid flowing be- 
tween two surfaces with a varying distance between them, 
see Fig. [1} will develop a heUcoidal component of veloc- 
ity with non-vanishing vorticity along the flow direction. 
This can be checked explicitly for the exactly solvable 
flow in a converging channel (§ 23 of Ref. 1). 

Another consequence of Eq. ([S]) is a possibility of 
separation of particles of different chirality in hydrody- 
namic flows. It has been observed in numerical simula- 
tions [ini HH mi HSl and recent experiments [15 . We 
note that according to Eq. ([9| in a stationary flow and in 
the linear approximation in the shear rate diVj, we have 
V^Wi — 0, and the flrst term in Eq. ([8| vanishes. Thus 
separation chiral isomers in the absence of external forces 
acting on the particles is possible either in non-linear or 
in non-stationary flows. 

In the practically important case of a stationary Cou- 
ette flow, the first term in Eq. (|8| vanishes for arbitrary 
Reynolds numbers, and the chiral current arises only due 
to orientation of the particles. The latter increases with 



the rotational Peclet number and saturates at Pe ^ 1. 
In this regime the chiral current becomes linear in the 
flow vorticity ui, Eq. (12). The linear dependence of jch 



on uj and saturation of the proportionality coefficient at 
Pe ^ oo has been observed numerically in Refs. [TOJ [T^ 

Separation of particles by chirality can also be achieved 
by subjecting the particles to an external circularly po- 
larized electric or magnetic field. The orientation of the 
particles along the field (e.g. due to the presence of a per- 
manent electric or magnetic dipole moment or anisotropy 
of the polarization matrix) will cause their rotation rel- 
ative to the surrounding fiuid. This will produce a chi- 
ral flux along the circular polarization axis due to the 
propeller effect [2l4|22j (similarly, a stationary electric or 
magnetic field will cause separation of particles by chiral- 
ity in a rotating fluid) . We also note that there are other 
mechanisms of chiral current generation which do not in- 
volve transfer of angular or linear momentum from the 
ac-field to the particles [23]. Chiral separation by circu- 
larly polarized magnetic field has been recently observed 
in the experiments of Refs. [53] The full quantitative 
analysis of this effect is beyond the scope of this work. 
Here we restrict the treatment to the experimentally rel- 
evant regime of strong and slowly varying fields, where 
thermal fluctuations can be neglected and the particles 
are fully polarized along the instantaneous electric field. 
In this case the problem is of purely geometric nature. 
The chiral current becomes independent of the viscosity 
of the fluid and can be expressed in terms of the Berry 
adiabatic connection [20]. Below we express this adia- 
batic connection in terms of the resistance matrix of the 
particle The latter relates the external force F and 
torque r exerted on the particle to the linear velocity 5-v 
and angular velocity Su) relative to the fluid. 



K C\f6w 
C h )\5ui 



(13) 



Here we chose the origin of the reference frame at the 
reaction center, so that the coupling tensor C is sym- 
metric 18J. Since a uniform electric field exerts no force 
on the particle we immediately obtain the relations. 



5w = -li-^Cduj, 



(14a) 
(14b) 



where we introduced the notation 17 = f2 — CK~^C . 

The orientation of the particle relative to the lab frame 
is described by the three Euler angles, (j), 9 and ijj [26]. We 
choose the axes of the body reference frame, Xi, X2, x^, so 
that cca points along the dipole moment of the particle. 
For fully polarized particles the latter points along the in- 
stantaneous direction of the electric field. Thus 9 and 
coincide with the polar angles of the electric field vector. 
The value of ^ remains undetermined because the parti- 
cle can be rotated by an arbitrary angle about without 
changing its polarization energy. When the orientation 
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of the electric field changes with time the particle ori- 
entation angle about the instantaneous field direction, 
tp{t), also changes. Its time evolution can be determined 
from the condition that projection of the torque onto the 
X3 axis must vanish. This is clear because the torque 
-r = d X E is perpendicular to the dipole moment d, 
which points along 2:3. Writing Eq. (14b I in the body 
frame, fisio;,; — 0, and expressing LUi in terms the Eu- 
ler angles, Wi = sin sin -0 + 9 cos ip, etc. we obtain 
dil; = Aj,d(l) + Agd9 where Aj, and Ag play the role of the 



adiabatic connection components and are given by 

A^ — — cos 9 — ^^i^ (^31 sin i/j + ^1^2 cos i/") , (15a) 
^33 ^ ' 



(15b) 



Ag — ~ - — ( r^si COS i\) + sin i/' 

^3\ 



This defines '0 in terms of 9 and The displacement of 
the particle can be obtained from Eq. ( 14a I . Instead of 



presenting the general formulae we focus on the practi- 
cally relevant case of an ac-field of frequency loq circularly 
polarized in the xy plane: 9 = 7r/2, (j) = uJot. In this case 
it is easy to see that the average velocity along the x 
and y axes vanishes. For the average ^-component of the 
velocity an elementary calculation gives, 



= xR^o 



Wo 



2n 



Tr 



33 



-^33 











-^33 





^31 


^32 






(16) 

where the resistance tensor is expressed in the body 
frame, and x can be viewed as the dimensionless mea- 
sure of the particle chirality. Since the coupling tensor C 
changes sign under inversion it is clear that particles of 
opposite chirality will move in opposite directions along 
the z axis. By order of magnitude the chiral separation 
velocity is v^^ ^ Ri^o, which is consistent with the recent 
experimental findings [241 ES] . 

In the regime where the polarization energy in the elec- 
tric field is smaller than the temperature the chiral cur- 
rent is reduced compared to the above estimate. The 
leading contribution at weak fields is proportional to the 
intensity of the ac-radiation [25] . 

So far we discussed the case where the suspended par- 
ticles consist of the opposite enantiomers of a single 
species. However, the effects considered above exist even 
in suspensions of particles of completely random shape in 
a non-chiral liquid. In this case the definition of chiral- 
ity requires clarification. For example one can define the 
chirality of a particle by considering the direction of its 
motion in a hydrodynamic shear flow or under the action 
of an ac-electromagnetic field. Thus the same individ- 
ual particle can exhibit different chirality with respect to 
different external perturbations. 

The set of Eqs. still holds for a suspension of 

random particles. In this case the auxiliary quantities 



n'^^ , and j '^^ are defined in terms of the chiral component 
of the stress tensor Eq. ^ and correspond to quantities 
averaged over the random shape of the particles. 

Finally, we note that symmetry allows contributions to 
jc/i that are proportional to the external magnetic field B, 
for example j'' cx n (yrfB. We believe that these ef- 
fects are of fluctuational origin similar to those discussed 
in Refs. IH |5l El [7] and do not study them in this work. 
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